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Abstract. We have developed a high-precision code which solves the kinematic dynamo problem both for given rotation law 
and meridional flow in the case of a low eddy diffusivity of the order of 10 11 cm 2 /s known from the sunspot decay. All our 
models work with an a-effect which is positive (negative) in the northern (southern) hemisphere. It is concentrated in radial 
layers located either at the top or at the bottom of the convection zone. We have also considered an a-effect uniformly distributed 
in all the convection zone. In the present paper the main attention is focused on i) the parity of the solution, ii) the form of the 
butterfly diagram and iii) the phase relation of the resulting field components. If the helioseismologically derived internal solar 
rotation law is considered, a model without meridional flow of high magnetic Reynolds number (corresponding to low eddy 
diffusivity) fails in all the three issues in comparison with the observations. However, a meridional flow with equatorial drift 
at the bottom of the convection zone of few meters by second can indeed enforce the equatorward migration of the toroidal 
magnetic field belts similar to the observed butterfly diagram but, the solution has only a dipolar parity if the (positive) a-effect 
is located at the base of the convection zone rather than at the top. We can, therefore, confirm the main results of a similar study 
by Dikpati & Gilman (2001). 
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1. Introduction 

In an early paper Roberts ( 1 972) discussed the excitation condi- 
tions of distributed shell dynamos. The a-effect, antisymmetric 
with respect to the equator, was considered in an outer shell 
(xi < x < 1, x = t/Rq), the dynamo was embedded in 
vacuum and the differential rotation dfl jdx was considered as 
uniform throughout the outer shell. He found for positive dy- 
namo numbers (a nort h ■ dSl/dx > 0) that dipoles are more eas- 
ily excited than quadrupoles for thin convection zones whilst 
for deeper zones the quadrupoles are more easily excited (see 
Fig-[|. 

These models are rather rough. The differences between the 
dipolar solutions and the quadrupolar solutions always happen 
to be small, so that the important question concerning the lati- 
tudinal symmetry, that is the parity, did not seem to be a serious 
problem. 

Roberts & Stix (1972) computed a 2 i?-dynamos with a ro- 
tation law similar to the solar rotation law known nowadays by 
helioseismology. The results are given in their Fig. 4b. Without 
latitudinal differential rotation (J3 = 0) but with dfi/dx\ e ^ > 
the solution with quadrupolar (symmetric) parity has a lower 
eigenvalue than the solution with the dipolar symmetry (oppo- 
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site to the Steenbeck-Krause (1969) model with d!l/dx\ c ^ < 
0, see Fig. 3 of Roberts & Stix). The inclusion of the latitudinal 
shear produces still higher dynamo numbers but the difference 
between quadrupolar parity and dipolar parity grows. 

Kohler (1973) thus considered only the excitation of modes 
with prescribed dipolar parity. For positive a-effect in the 
northern hemisphere and with an outwards increasing i7, only 
a poleward drift of the toroidal magnetic field belts was found, 
contrary to the observations. 

Moss & Brooke (2000), in order to produce the observed 
equatorward migration of the toroidal fields adopt the solar ro- 
tation law and a negative northern a-effect in the bulk of the 
convection zone. The dipolar solutions are only slightly eas- 
ier to excite than the quadrupolar ones (dipole: C a — —3.20, 
quadrupoleC Q = —3.25). The parity problem does not seem to 
exist for negative dynamo numbers (see also Fig. 3 in Roberts 
& Stix 1972). The situation, however, completely changes if a 
positive northern a-effect is considered (Roberts & Stix 1972, 
Fig. 4; Moss 1999). In this case one has to deal with a parity 
problem for the solar dynamo. 

Following Dikpati & Gilman (2001) we shall here study 
the parity problem for solar dynamos in particular for dynamos 
with rather small eddy diffusivities so that the meridional flow 
plays an important role in advecting the toroidal magnetic field 
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Fig. 2. The stream function ip of the meridional circulation used 
in the model computations. It must be negative if the circulation 
shall proceed equatorwards at the surface. 



Fig. 1. The ratio of the dynamo numbers for equatorially an- 
tisymmetric (dipoles) to equatorially symmetric (quadrupoles) 
solutions for the shell dynamo of Roberts (1972). At x\ = 0.65 
both values are equal. For thinner shells the most efficient dy- 
namo is of dipolar type while for ticker shells is of quadrupolar 
type- 



belts (Wang et al. 1991; Choudhuri et al. 1995). As the rotation 
law can be considered as given (known from helioseismology), 
we are free to vary the location of the a-effect, so that models 
are assumed to have a positive a-effect both at the top and at the 
bottom of the convection zone, as well as in the full convection 
zone. 

The inclusion of the meridional flow u m has a strong im- 
pact on the mean-field dynamo when the eddy diffusivity 77^ 
is low. In particular, for ?7t= 10 n cm 2 /s , as is known from 
the sunspot decay, the magnetic Reynolds number Rm = 
u m R0/r]T reaches values of the order of 10 3 for a flow of 
10 m/s. As a consequence, depending on the localization of 
the dynamo wave, a dramatic modification of both magnetic 
field configurations and cycle period is expected. This possibil- 
ity has recently been a subject of intense numerical investiga- 
tion (Dikpati & Charbonneau 1999; Kiiker et al. 2001), where 
it has been shown that solutions with high magnetic Reynolds 
number provide correct cycle period, butterfly diagrams, mag- 
netic phase relations and sign of current helicity, by means of a 
positive a-effect in the north hemisphere. 

In this investigation we study how the presence of the flow 
and the location of the turbulent layer affect the parity mode 
selection and the cycle period. In this respect we show that, for 
high magnetic Reynolds numbers of the flow, quadrupolar field 
configurations are more easily excited than the dipolar ones if 
there is no a-effect below r/R Q m 0.8. 



2. Basic equations 

In the following the dynamo equations are given with the inclu- 
sion of the induction by meridional circulation. Axisymmetry 
implies that the mean flow in spherical coordinates is given by 



In our formalism the magnetic field reads 
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where A is the poloidal-field potential and B is the toroidal 
field. The dynamo equation may be written in the form 

dB 

— = rot(tt x B + aB) - rot (rot y/rf^B)) , (3) 

which includes the diamagnetism due to non-uniform turbu- 
lence. In fact it has been shown (Kitchatinov& Riidiger, 1992) 
that the diffusive part of the mean turbulent electromotive force 
reads — i]^TotB— Vr/TxB/2 = — ^/rfrTot^rjrB. The second 
term is due to the turbulent diamagnetism. If there are strong 
gradients of turbulence intensity, this term will dominate the 
transport of the mean magnetic fields. 

As usual, the meridional circulation is derived from a 
stream function, so that the condition div pit = is automati- 
cally fulfilled. Then a series expansion in Legendre polynomi- 
als is introduced, as described in Riidiger (1989): 
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A one-cell meridional circulation is described by 
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where ip is the usual stream function. A positive ip describes a 
cell circulating clockwise in the northern hemisphere, i.e. the 
flow is polewards at the bottom of the convection zone and 
equatorwards at the surface. For a negative ip the flow is, as 
is observed, polewards at the surface. In order to keep the flow 
inside the convection zone, the function ip must be zero at the 
surface and at the bottom of the convection zone. In order to de- 
fine the strength of the flow, we shall use the values u m of the 
meridional circulation at the bottom of the convection zone. 
The rotation law for the solar dynamo can be considered as 
given by the helioseismic observations, in particular we have 
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used the analytical model of Dikpati & Charbonneau (1999) 
which is characterized by the existence of a steep subrotation 
profile in the polar region with a thickness of about 0.05 solar 
radii (see Markiel & Thomas 1999). We then define the profile 
of as 

VT = Vc + \{m - tfc)(l + erf(40(x - 0.7))), (8) 

where x — r/R Q is the fractional radius, erf denotes the error 
function, r\ t is the eddy diffusivity, and T) c the magnetic diffu- 
sivity beneath the convection zone. The factor 40 defines the 
thickness of the transition region to be 0.05i?©, and we have 
adopted both r) c /rj t = 0.1 and rj c /rj t = 0.02. 

There is an increasing evidence for an a-effect in the con- 
vection zone rather than in the overshoot region. Observations 
of the turbulent current helicity by Seehafer (1990), Pevtsov 
etal. (1995), Abramenko etal. (1996) as well as Bao & Zhang 
(1998) always lead to negative values in the northern hemi- 
sphere and positive values in the southern hemisphere. If the 
observations indeed concern the current helicity of the fluctu- 
ations then after Radler & Seehafer (1990) the azimuthal a- 
effect should be positive. Also Brandenburg (1999) with new 
numerical simulations found (highly noisy) positive correla- 
tions between the turbulent EMF and the mean magnetic field, 
i.e. positive a-effect for the northern hemisphere. Our SOCA 
theory of Lorentz force-driven turbulence in a stratified rotating 
plasma also leads to negative current helicity and positive a- 
effect for the northern hemisphere (Riidiger et al. 2001). Lastly, 
with their high-detailed simulation, Ossendrijver et al. (2001) 
show that indeed far beyond the limits of the SOCA theory the 
results are confirmed: in the northern part of the solar convec- 
tion zone the horizontal a-effect proves to be positive at vari- 
ance with the vertical a-effect and the current helicity. 

The a-effect is always antisymmetric with respect to the 
equator, so that we write 

a = a cos + erf(^^)) (l - erf(^-^))/4, (9) 

where ceo is the amplitude of the a-effect and x a , xp and d 
define the location and the thickness of the turbulent layer, re- 
spectively. Differently from the overshoot dynamo ceo is not 
assumed to change its sign in the bulk of the convection zone 
or in the overshoot layer. In Fig. |3| the eddy diffusivity profile 
and the a profile are shown for x a — 0.9, xr = 1 and d = 0.05 
which correspond to an a-effect located at the top of the con- 
vection zone. 

3. Mathematics and numerics 

In order to investigate the properties of the solar dynamo with 
such a low eddy diffusivity we solve the linear dynamo Eq. (|3|) 
with a finite-difference scheme for the radial dependence and a 
polynomial expansion for the angular dependence. In particu- 
lar, we have used the following expansions for the field: 

A(x,6) =e xt 5> n (z) pW(cos0), (io) 
B(x,9) =c xt ^b m {x) P£\eo S e), (11) 
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Fig. 3. The a-effect normalized function in the bulk of the con- 
vection zone (solid line) and the profile of the eddy diffusivity 
(dashed line). 
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Fig. 4. Truncation dependence of a typical high-Reynolds num- 
ber solution. The quantities ma3i x {b ri } and max x {a n } as a 
function of n are shown, in the left and in the right panel, re- 
spectively. 

where A is the (complex) eigenvalue, n = 1, 3, 5, . . . and m = 
2,4,6,... for antisymmetric modes, and n ^ m for symmetric 
modes. Vacuum boundary conditions at the surface are then 
translated into 

-^ + (n + l)a n = b m = 0. (12) 
dx 

In the interior at x = x\ — 0.6 we have instead set 

db rn 

x— \-b m = a n = 0, (13) 

dx 

which imply perfect conductor boundary conditions for the 
field. 

By substituting (|l0|) and ( pi] ) in (]3|) one obtains a set of in- 
finite o.d.e. that can be conveniently truncated in n when the 
desired accuracy is achieved. The system is in fact solved by 
means of a second order accuracy finite difference scheme and 
the basic computational task is thus to numerically compute 
eigenvalues and eigenvectors of a block-band diagonal real ma- 
trix of dimension M x n, Al being the number of mesh points 
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and n the number of harmonics. This basic algorithm is embed- 
ded in a bisection procedure in order to determine the critical 
a-value needed to find a purely oscillatory solution. This value 
is accepted when the dimensionless quantity 5fte(A)i?Q/?7T is 
no greater than 1CP 3 . 

We have first tested the accuracy of our code in simple cases 
where the eigenvalues and eigenfunctions are known (decay 
modes with a — and simple analytic solutions with con- 
stant a). Then, we have checked most of the solutions dis- 
cussed in the literature (no flow) and we found good agree- 
ment. In fact we have rapid convergence in n for simple (no- 
flow) a/2-dynamo, and we could confirm most of the findings 
by Steenbeck & Krause (1969) and Roberts (1972). When the 
flow is present, the number of harmonics needed in order to 
get convergence generally increases with the Reynolds number 
and we decided to truncate the series (|l2[|lj) when the maxi- 
mum value over x of the n-th harmonic drops by roughly three 
orders of magnitude, as shown in Fig. |j. 

4. Numerical results 

We performed an extensive numerical exploration of the pa- 
rameter space, allowing simultaneous variations of x a , r] c /r)t, 
u m and angular dependence of the a-effect, considering both 
cos and cos sin 2 9 functions in In all models discussed 
in the following the rotation law is given by helioseismology 
and the eddy diffusivity is fixed to be ?/t = 10 11 cm 2 /s. The 
drift amplitude is considered as a free parameter of the order 
of 1 m/s up to 10 m/s. In some cases we have investigated the 
field configuration for very strong flow up to 20 m/s. 

4. 1. Alpha-effect in the entire convection zone 

We start with a model where the (positive) a-effect exists 
through the whole convection zone (x a = 0.7, xp — 1, 
d = 0.05). In Table |l]the results are given. Suffix A denotes the 
(dipolar) solutions with antisymmetry with respect to the equa- 
tor and suffix S denotes the (quadrupolar) solutions with sym- 
metry with respect to the equator. The drift amplitude varies 
between 2 cm/s and 10 cm/s. For small flow the dipole-solution 
occurs with the lowest a-effect amplitude, but for higher values 
of the flow the quadrupole solution is more easily excited. The 
oscillation period of the linear dynamo is also given: for high 
value of the flow quadrupole solutions have cycles which are 
longer than the corresponding dipole ones. This is a property 
that we have verified for all cases in the rest of paper, namely 
that in the high Reynolds number regime, lowest a-value solu- 
tions have longer periods than the corresponding opposite par- 
ity solutions. 

It should also be noticed that, with the small value of eddy 
diffusivity we have chosen, the cycle period becomes rather 
long compared with the 22-years of the Sun. As we know from 
the theory of the overshoot dynamo, the inclusion of the nonlin- 
ear feedback of the magnetic field (via a-quenching) reduces 
the periods to more realistic values (Riidiger & Brandenburg 
1995). What we are interested in, is mainly the influence of the 
a-effect profile and the drift amplitude on the resulting parity 
of the solutions. Figures @ and H show the magnetic topology of 



the dynamo with 3 m/s drift amplitude. We find (for the lowest 
eigenvalue) a solution with antisymmetry with respect to the 
equator. The toroidal field belts are concentrated at the bottom 
of the convection zone, and the poloidal field exhibits, close to 
the surface, a rather small-scaled structure (Fig. |J). The toroidal 
field belts migrates equatorwards but, however, the maximal 
field strength occurs in the polar region. 




Fig. 5. Alpha-effect in the entire convection zone: Toroidal 
(left) and poloidal (right) antisymmetric (dipolar) field con- 
figuration for u m = 6 m/s at t = 0. The a-effect is located 
between x a = 0.7 and xp = 1. The solution is plotted as B^- 
isocontours and poloidal magnetic field lines. Solid contours 
(field lines) correspond to positive toroidal fields (pointing into 
the plane of the paper) and dashed contours to negative toroidal 
field (pointing out of the paper) . 



Table 1. Critical a-values (cm/s) and periods P (yrs) for 
models with a-effect in the entire convection zone. Bold is 
used when the dipolar solution has the smallest a-value. The 
quadrupolar symmetry is denoted by "S" ("symmetric") and 
the dipolar parity symmetry by "A" (antisymmetric). 
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4.2. Alpha-effect at the top 

We now consider models where the a-effect is located at the 
top of the convection zone as discussed in the models of Kiiker 
et al. (2001). The results of the simulations confirm the basic 
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Fig. 6. Alpha-effect in the entire convection zone: Butterfly di- 
agram (top) and field phase relation (bottom) of a dynamo with 
critical turbulence ao = 2.46 cm/s, u m = 6 m/s. Black colors 
means negative B r B^. The a-effect is located between x = 0.7 
and x = 1. 
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Fig. 7. Alpha-effect at the top: Butterfly diagram (top) and field 
phase relation (bottom) of a dynamo with critical turbulence 
ao = 10.8 cm/s, u m = 10 m/s, x a — 0.9, xp = 1. Black 
colors means negative B r B<p. 

features of the advection-dominated dynamo, namely that, for 
a flow of few m/s and a low diffusivity, the butterfly diagram 
shows the correct equatorward migration of the toroidal field, 
and the phase relation of the magnetic fields is mostly negative 
as shown in Fig. ^ We found that, as far as the parity model 
selection is concerned, variations of the ratio r] c /rjt are not par- 
ticularly significant, while a functional dependence of the type 



cos 9 sin 2 8 for the a-effect disfavours the appearance of dipo- 
lar field configurations. We have then considered variations of 
the thickness of the a-layer. The results of this investigation 
are summarized in Fig. ^| and in Tables || and ||. For a slow 



Table 2. Alpha-effect at the top (thick layer): Critical a-values 
(cm/s) and periods (yrs) for various values of the flow (m/s), for 
both antisymmetric (A) and symmetric field (S) configurations 
with x a = 0.8, xp = 1, rjc/vt = 0.02. 
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Table 3. Alpha-effect at the top (thin layer): Critical a-values 
(cm/s) and periods (yrs) for various values of the flow (m/s), 
for both antisymmetric and symmetric field configurations with 
x a — 0.9, xp — 1, and rj c /r]t = 0.02 . The dipole solutions 
with lowest critical a- value are shown in bold. 
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flow, both quadrupoles and dipoles have similar excitation con- 
ditions. For intermediate values of the flow amplitude, it is pos- 
sible to have a different dynamo mechanism working and the 
solutions are stationary. In particular, as it is shown in Tables || 
and ||, dipole solutions may have a smaller critical ao _vaiue m 
this case. 

In the high-Reynolds numbers regime, quadrupole fields 
are more easily excited and the cycle periods drastically reduce. 
We confirm the findings by Dikpati & Charbonneau (1999) as 
far as the global structure of the field and the parity problem 
(Dikpati & Gilman 2001) is concerned. For a flow of 10 m/s 
the field geometry is shown in Fig. ^ and the mode with the 
lowest eigenvalue is a quadrupole. 
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Fig. 8. Alpha-effect at the top: The ratio cxa/cxs of the critical 
a-numbers for antisymmetric and symmetric field configura- 
tions is shown (as in Fig. [l]) as a function of the flow speed u m 
in the panels above. The period (continuous line for the anti- 
symmetric modes, dashed line for symmetric modes) is shown 
in panels below. The presence of vertical lines indicates that the 
solution is stationary 



4.3. Alpha-effect at the bottom 

This picture drastically changes if the a-effect operates at the 
bottom of the convection zone. We have considered a thin and 
a thick a-layer. In the first class of models x a = 0.7 and xp = 
0.8 while in a second class x a = 0.7 and xp = 0.9. 

Figures |To| and |ll] show an example of a model of the first 
class. Again, the toroidal field is concentrated at the bottom of 
the convection zone, where, however, the highest field ampli- 
tudes occur in the polar regions. The diagram for B r ■ shows 
the dominance of the negative sign. 

In Fig. |l2] parity and cycle periods of the solution for the 
model with a thin a-layer at the bottom of the convection zone 
are given. The solution with the dipolar symmetry has the low- 
est a-value and is the stable one. There is no heavy parity prob- 
lem if the a-layer is located at the base of the convection zone 
(see Dikpati & Gilman 2001). However, we can also notice 
from Fig. [1^ (lower panel) that in this case the oscillation pe- 
riod of the dipoles is longer than for the quadrupoles, while, 
for slow flows, the cyclic behaviour of the dynamo solution 
disappears, so that only meridional flows with amplitudes ex- 
ceeding 3 m/s are here relevant. In this case the dipole solutions 
do not match the 22-year cycle period of the Sun (although the 





Fig. 9. Alpha-effect at the top: Toroidal (left) and poloidal 
(right) symmetric (quadrupolar) field configuration for ao = 
8.95 cm/s, x a = 0.9, xp = l,d = 0.02, 17,5/77* = 0.1, u m = 10 
m/s at t = 





Fig. 10. Alpha-effect at the bottom: Toroidal (left) and poloidal 
(right) antisymmetric (dipolar) field configuration for a= 0.95 
cm/s, u m = 6 m/s at t = 0. The a-effect is located between 
x a — 0.7 and xp = 0.8. 



quadrupoles do). The overall result is that the dipolar solutions 
are always more easily excited and also the butterfly diagram 
shows the correct characteristics. 

In the second class of models the dynamo mechanism was 
of the same type, as we could infer from the field configuration 
and from the butterfly diagram (not shown). However in this 
case the quadrupolar solution was always the most easily ex- 
cited. From these results we can deduce that the region where 
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Fig. 11. The same as in Fig. [7] but for the model presented in 
Fig. [k]. The cycle period is 120 years. 

the a-effect produces more easily dipolar field configuration, 
in the advection-dominated regime, is below x s» 0.8. 

4.4. Very thin a -layer 

In this case we have set x a = 0.7, xp — 0.75, d — 0.02 and the 
relevant results are summarized in Table ^j, where it is possible 
to notice that for a thin a-layer at the bottom, the dipolar field 
configuration is the most easily excited. We noticed that if no 
flow is present there is no migration towards lower latitudes of 
the magnetic field. However for moderate flow the maximum 
of the toroidal field is in the equatorial region, as is possible to 
see in Fig.|l3|and Fig.[l4]. For stronger flow the period matches 
the observed one as shown in Table [|. In this case dipolar and 
quadrupolar solutions have greater dynamo numbers than in the 
slow flow case, but the dipolar field configuration is always the 
favoured one. The field configuration and the butterfly diagram 
are shown in Fig. Conversely, a very thin a-layer lo- 

cated at the top of the convection zone produces a much larger 
dynamo number and the solution is always symmetric as shown 
in Table @ 

5. Conclusions 

A kinematic solar dynamo has been studied with the flow pat- 
tern taken from observations (rotation law) and theory (merid- 
ional flow). An eddy diffusivity of 10 11 cm 2 /s provides us with 
a value consistent with the sunspot decay. With such a small 
value the magnetic Reynolds number for a meridional flow of 
(say) 10 m/s reaches values of order of 10 3 , so that the dynamo 
can really be called advection-dominated. The meridional flow 
at the bottom of the convection zone is supposed to drift equa- 
torwards (and polewards at its top) as can be deduced by the 
mean-field theory of the differential rotation in an outer convec- 
tion zone (Kitchatinov & Rudiger 1999; Miesch et al. 2000). 
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Fig. 12. The same as in Fig. ||but for models with the a-effect 
located at the bottom of the convection with x a = 0.7 and 
xp — 0.8. Note the clear dominance of dipolar solutions indi- 
cated by their much smaller dynamo numbers. 

Table 4. Critical a-values (cm/s) and periods (yrs) for a very 
thin a-layer located between x a = 0.7 and xp — 0.75 for 
various values of the flow (m/s), for both antisymmetric and 
symmetric field configurations. In the last line the same a-layer 
is located between x a = 0.95 and xp = 1 and the solution 
is clearly of the symmetric type. Also notice the much larger 
dynamo number in this latter case. 
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When this happens, the meridional flow advects the field equa- 
torwards producing a butterfly diagram of the observed type, 
which would not occur with i) a positive a-effect (in the north- 
ern emisphere), ii) the standard rotation law known from helio- 
seismology and iii) no meridional flow (Choudhuri et al. 1995; 
Dikpati & Charbonneau 1999; Kiiker et al. 2001). As Dikpati 
& Gilman (2001) have stressed, these models may encounter 
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Fig. 13. Thin Alpha-effect at the bottom: Toroidal (left) and Fig. 15. Alpha-effect at the bottom for a very thin layer: 
poloidal (right) antisymmetric (dipolar) field configuration for toroidal (left) and poloidal (right) antisymmetric (dipolar) field 
a= 4.36 cm/s, u m = 3 m/s at t = and x a = 0.7, xp = 0.75. configuration for a= 16.61 cm/s, u m = 12 m/s at t = and 

x a — 0.7, Xf3 — 0.75 . 




phase angle 




phase angle 



Fig. 14. The same as in Fig. ^ but for the model presented in 
Fig. fijl. The cycle period is 74 years. 




4 6 
phase angle 




Fig. 16. The same as in Fig. [7] but for the model presented in 
Fig. O. The cycle period is 23 years. 



a problem with the parity of the solution since for most of the 
models the quadrupolar solution is the stable one. We have used 
this striking fact to discuss the effect of the radial distribution 
of the a-effect on the parity model selection. Our conclusion 
is that a thin a-layer located below r/R Q » 0.8 selects mod- 
els with correct parity property, butterfly diagram phase rela- 
tion and cycle periods close to the observed one. This confirms 
the findings of Dikpati & Gilman (2001), providing further ev- 
idence for a tachocline a-effect as a promising candidate for 
understanding the dynamo mechanism operating in the Sun. 
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